;***************************************************************
;discription：初级生产力-NPP计算
;
;param1: nrows(input) 行
;param2: ncols(input) 列
;param1: apar(input)植物的光合有效辐射 单位为gc*m^-2*month^-1
;param2: e(input) 实际光能利用率 单位为gc*mj^-1
;param3: npp(output) 区域植被初级生产力
;
;author: wangqiang( wensten )  2015-09-15  chengdu
;***************************************************************
PRO PRO_CASA_NPP, nrows, ncols, apar, e, npp=npp
  ;
  FOR yindex=0,nrows-1 DO BEGIN
    FOR xindex=0,ncols-1 DO BEGIN
      ;获取光合有效辐射和实际光能利用率
      apar_cellvalue=apar[yindex,xindex]
      e_cellvalue=e[yindex,xindex]
      
      IF (apar_cellvalue EQ -9999.0)||(e_cellvalue EQ -9999.0) THEN BEGIN
        ;无数据区
        npp[yindex,xindex] = -9999.0
      ENDIF ELSE BEGIN
        res = apar_cellvalue * e_cellvalue
        ;去除噪音
        res= res > 0
        npp[yindex,xindex]=res
      ENDELSE
      
    ENDFOR
  ENDFOR
  
END